Compton Heated Outflow from CDAFs 
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ABSTRACT 

Convection-dominated accretion flows (CDAF) are expected to have a shal- 
lower density profile and a higher radiation efficiency as compared to advection- 
dominated accretion flows (ADAF). Both solutions have been developed to ac- 
count for the observed properties of the low luminosity, high temperature X- 
ray sources believed to involve accretion onto massive black holes. Self-similar 
CDAFs also have steeper poloidal density gradients and temperatures close to 
the virial temperature. All these characteristics make CDAFs more capable 
of producing polar outflows driven by Compton heating as compared to other 
classical accretion disks. We investigate the conditions for producing such out- 
flows in CDAFs and look for the mass accretion rate, or, equally, the luminosity 
of CDAFs for which such outflows will exist. When the electron temperature 
saturates around 10 11 K at the inner region, polar outflows are probable for 
8 x 10~ 7 < L/L E < 4 x 1CT 5 , where L E is the Eddington luminosity. Out- 
flows are well collimated with small opening angles. The luminosity range for 
which outflow solutions exist is narrower for lower electron temperature flows and 
disappears completely for electron temperature <6x 10 9 K. When the magnetic 
field is present, we find that outflows are possible if the magnetic field is less than 
from 10% to 1% of the equipartition field. We also find that outflows are more 
likely to be produced when the viscosity parameter a is small. The tendency 
for jet-like collimated outflows for these solutions is presumably astrophysically 
relevant given the high frequency of jets from AGNs. 
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Subject headings: accretion, accretion disks — black hole physics — quasars: gen- 
eral — X-rays: general 

1. Introduction 

Advect ion-dominated accretion flows (ADAF) nicely complement the classic thin disk 
accretion flows (Shakura & Sunyaev 1973; Ichimaru 1977; Rees et al. 1982; Narayan & Yi 
1994, 1995; Abramowicz et al. 1995), and have been successfully applied to variety of cosmic 
objects, from galactic X-ray binaries to diffuse X-ray background (see Narayan, Mahadevan, 
& Quataert 1999 for review). However, analytic studies of ADAFs indicated the convectively 
unstable nature of ADAFs (Narayan & Yi 1994, 1995a; see Begelman & Meier 1982 for 
radiation-dominated ADAF), which has subsequently been proved in a series of numerical 
simulations (Igumenshchev, Chen, & Abramowicz 1996; Igumenshchev & Abramowicz 1999, 
2000; Stone, Pringle, & Begelman 1999; Igumenshchev, Abramowicz, & Narayan 2000). 
Especially, the numerical studies of ADAFs by Igumenshchev & Abramowicz (1999, 2000) 
show that an ADAF becomes a convection-dominated accretion flow (CDAF) whenever the 
viscosity parameter a < 0.1. Further analyses of self-similar CDAF solutions show their 
unique properties (Narayan, Igumenshchev & Abramowicz 2000, hereafter NIA; Quataert 
& Gruzinov 2000, hereafter QG): the density varies as p oc i?" 1 / 2 (R is the radius), much 
flatter than the usual i?™ 3 / 2 in an ADAF or in spherical accretion. Correspondingly, the 
mean radial velocity varies as v oc R~ 3 / 2 , compared to R~ x l 2 in ADAF or in spherical flow. 
Energy generated at the inner part of the flow is transported to the outer part by convection. 
A CDAF has perhaps as much resemblance to the rotating stellar envelope of a massive star 
as to the usual accretion flow. 

The other aspect of the two-dimensional nature of ADAF or CDAF solutions, often 
neglected, is the interaction between the outgoing radiation produced at smaller radii with 
the inflowing gas in the outer part of the flow. In optically thick stars this interaction plays 
a vital role in establishing the equilibrium states. The radiative interaction also plays a very 
important role in pure spherical accretion flows (Ostriker et al. 1976; Cowie, Ostriker, & 
Stark 1978; Wandel, Yahil, & Milgrom 1984; Park 1990a, 1990b; Nobili, Turolla, & Zampieri 
1991; Zampieri, Miller, & Turolla 1996; Ciotti & Ostriker 1997, 2001). Park & Ostriker 
(1999, 2001) studied the same radiative interaction in the context of the ADAF solution 
and found that a polar outflow can be generated through Compton heating of electrons by 
high-energy photons emitted by the inner, hot part of the flow. The winds generated by the 
processes in the papers listed above are not momentum driven. Rather, they are caused by 
overheating of the gas in the slowly moving, low density polar regions. Outflows may also be 
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generated from ADAF by other hydrodynamic processes (Narayan & Yi 1995; Xu & Chen 
1997; Blandford & Begelman 1999). 

In this work, we study the conditions for CDAFs to develop radiation driven outflows. 
We adopt the self-similar CDAF solution as the background flow structure (NIA; QG). The 
treatment in this paper is two dimensional, adopting the angular profile of NIA and QG, 
except that the radiation field is simplified and treated as spherically symmetric. 



2. Flow Properties 

2.1. Density and Temperature 

Multi-dimensional numerical simulations and analytic analyses show that density and 
temperature profiles of CDAF follow a self-similar form in radius (NIA; QG; Ball, Narayan 
& Quataert 2001, hereafter BNQ): 

p(r) oc r-- 1/2 (1) 
Ti(r) oc r-\ (2) 

where r is the radius in units of the Schwarzschild radius, R s = 2GM/c 2 = 3.0 x 10 5 mcm, 
and m is the black hole mass in solar units. This shallow radial dependence corresponds to 
marginal stability to convection. QG further showed that the density should follow a power 
law in sin-$ (■$ is the angle from the pole) while the temperature is constant on spherical 
shells. So we shall assume that the density is given by 

p(r,tf) = p r~ 1/2 sin 2 $ (3) 
= Pout(r/r out )-^ 2 sin 2 (4) 

The flow extends from r out down to r in — 1, with the outer radius (analogous to the stellar 
photosphere) fixed by energy balance considerations. 

The electron temperature T e in an accretion flow is determined by the balance between 
various cooling and heating processes. In general, detailed studies of hot accretion flows 
(Narayan & Yi 1995b; Narayan, Barret, & McClintock 1997; BNQ) show that the electron 
temperature is equal to the ion temperature for T e < 10 9 K, and then, at smaller radii, 
due to highly efficient relativistic bremsstrahlung and synchrotron emission, it flattens to 
somewhere between 10 9 K to 10 11 K, depending on the amount of direct viscous heating to 
electrons. Hence, in this paper, we approximate the electron temperature as 

T e (r,tf) = — for r>n (5) 
r 

= T\ for r < n, (6) 
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where To ~ 10 12 K, T\ is some constant between 10 9 K and 10 11 K, and r\ = Tq/T\. 

From equation (j3J), we find that the total Thomson optical depth along a given direction 

i? is 

r es ~ -mr o ; 4 sin tf, (7) 

where fn = M / 'Mg = Mc 2 /Le is the dimensionless mass accretion rate. Our analysis, which 
is to be applied to low accretion rate solutions, will assume that the flow is optically thin to 
both scattering and redistribution processes. 



2.2. Radiation Field and Radiation Temperature 

The inner part of the accretion flow generally produces higher energy photons, while 
the outer part produces lower energy photons. Hence the radiation field at a given position 
is determined by the contribution from inner radiating shells plus the outer shells. Although 
in each shell, the density, and therefore, the emissivity (i.e. cooling function) is a function 
of poloidal angle the energy density is more uniform over angle than is the emissivity. 
So, we will simplify the calculation of radiative transfer by assuming each radiating shell is 
homogeneous, using an angle average over the sphere. 

Therefore, the radiation energy density at given r, Ex{r), is given by 



cE x (r) = R s 



i r r r r out i rp I r 

— ^ / A(r')47ir' 2 dr' + / A(r')- In J— — dr' 
4vrr 2 L J r r \ r' - r 



where the first term is the contribution from inner radiating shells and the second term from 
outer radiating shells (see Appendix A). The emissivity per unit volume is denoted as the 
cooling function A, which will be described in §3.1. 

For relativistic un-Comptonized bremsstrahlung, with A oc p 2 T e and CDAF scaling p oc 
r -1 / 2 and T e oc r -1 , the incremental luminosity dLx(r)/dr = constant, and the luminosity 
at a given radius r has the largest contribution from the outermost shell. Nevertheless, 
equation gives a peak energy density in the inner region (cf. dashed line in Fig. 1) and 
a significant radiation pressure gradient. In addition, when there is strong Comptonization, 
the contribution from the inner region can be significantly enhanced. 

Moreover, the amount of Compton heating is determined by the luminosity times the 
photon energy, and is largely dominated by the inner hot region. The radiation temperature, 
defined as kTx(r) being the energy- weighted mean photon energy, at a given radius r is 
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determined by the equation 

cT x (r)E x (r) = R s 



1 pr pr out l I i _i_ 

- J T x (r>)A(r>)4nr"dr> + J T x (r')A(r')- In ^ -p-dr' 



Airr 2 

(9) 

This quantity is more concentrated than E x (r), especially in the presence of strong Comp- 
tonization [see dotted curve in Fig. 1 for T x (r)}. 

In presence of an overheated wind, some region along the pole in CDAF will not be able 
to accrete matter, while the remaining equatorial region normally accretes (see §3.3). In 
most of the relevant parameter space, the overheated 'funnel' is very narrow, and we assume 
(with regard to energy generation) that the self-similar CDAF flow is filling the whole space 
including the polar region, omitting the small correction due to the solid angle Qw (~ 0.2 
sr) occupied by the outflowing wind/jet. 

The resulting profiles of radiation moments and the radiation temperature show the 
aforementioned characteristics of CDAFs. The luminosity, 

L(r) =R% i A(r')4irr' 2 dr', (10) 

increases at large radius (solid curve in Fig. 1), while the radiation temperature is roughly 
flat or decreasing slowly outward after peaking at some intermediate radius (dotted curve in 
Fig. 1). 



2.3. Outer Boundary 

The mathematical self-similar CDAF has zero mass accretion rate (NIA), and this, 
of course, corresponds to zero luminosity. However, the accretion flow must lose energy 
by radiative cooling, and this energy loss is very likely provided by the convective energy 
transport (NIA, BNQ). We follow BNQ to assume that the total luminosity of ADAF is 
proportional to the mass accretion rate: 

L(r out ) = r] c e c Mc 2 , (11) 

where the convective efficiency e c , determined by numerical simulations to be in the range 
~ 1CT 2 — 1CT 3 and i] c is the fraction of convected energy to be radiated away by electrons. In 
subsequent examples, we nominally adopt e c i] c = 1CT 2 (BNQ). We also tried smaller values 
of e c i] c or a scaling of e c i] c with m (cf. §3.3), and the results were qualitatively the same. 
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We determine the outer boundary for a given mass accretion rate rh by finding the 
radius r out at which equation ffTTl) is satisfied. This gives a somewhat different relation be- 
tween the total luminosity and the mass accretion rate from that in BNQ because we have 
used Comptonized relativistic bremsstrahlung (and synchrotron) while BNQ employ only 
un-Comptonized non-relativistic bremsstrahlung. However, the general characteristic of the 
relation is the same: higher mass accretion flows have smaller outer boundary radii. This 
is due to the special radial profile of density in CDAF solutions. Because of the weak den- 
sity gradient with radius, most of the luminosity contribution is from the outer part of the 
flow, which makes CDAF flows quite distinct from other accretion flows including ADAFs. 
The dimensionless total luminosity, L{r out )/ L E , is directly proportional to the dimensionless 
mass accretion rate rh = M/Me by equation (fTTj) while the luminosity accumulated from 
bremsstrahlung is very roughly proportional to rh 2 r out (with a higher power in rh if Comp- 
tonization is important). Since these two luminosities must be equal, a higher rh flow implies 
a smaller r out (see solid lines in Fig. 2). Therefore lower mass accretion rate CDAF solutions 
have lower luminosity, yet are more extended, if we treat e c r] c as fixed, independent of rh. 



The main cooling mechanism we consider in this work is Comptonized bremsstrahlung 
and Comptonized synchrotron. We first focus on the Comptonized bremsstrahlung. This 
process is not only the source of electron cooling but also the source of the Compton heating 
radiation field. 

The Comptonized bremsstrahlung cooling rate per unit volume can be reasonably ap- 
proximated by (Svensson 1982; Stepney & Guilbert 1983) 



3. 



Outflow 



3.1. Cooling and Heating 




(12) 



where 




for 6 e < 1 



for 6 e > 1 



(13) 



ee 



4^(44 - 3n 2 )6 3 J 2 (l + l.W e + 6 2 e - 1.250 e 5 / 2 ) for 6 e < 1 
e pn(1.1230 e ) + 1.2746] for 6 e > 1, 



(14) 
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9 e = kT e /m e c 2 , n is the electron (ion) number density, ot the Thomson cross section, otf the 
fine structure constant, and m e the electron mass. In this work, we assume pure hydrogen 
gas. 

The emitted bremsstrahlung photons are upscattered by inverse Compton scattering. 
The amplification of photon energy by single scattering is A = 1+A9 e + 168 2 . The probability 
of single scattering is P = 1 — exp(— r es ) (see e.g., Dermer et al. 1991). Therefore the mean 
amplification factor by Compton scattering would be 

r] = 1- P + PA= 1 + P(A-l). (15) 

The fraction 1 — P of photons remain unscattered while the fraction P of photons are 
upscattered to A times their initial energy. This prescription is valid only for a single 
scattering Comptonization. Since most CDAF flow considered here has r es < 1, we use 
this amplification factor for Comptonized bremsstrahlung. Although Dermer et al. (1991) 
provides a handy formula for rj, applicable to diverse regimes, it can give incorrect values for 
A ^> 1 and r es < 1, the main parameter regime of CDAF. Although r es is a hard-to-define 
quantity in a complex flow, we simply adopt r es = n e aTr in this work for the purpose of 
estimating the Comptonization. 

When PA ^> 1, Comptonization becomes saturated, and all photons that are not ab- 
sorbed are upscattered to 3T e , that is 77 — > i] sat = 38 e /x where x = hv/m e c 2 (Dermer et al. 
1991). Then the fully saturated Comptonized bremsstrahlung emission is 

r-36» e 



a sat 
A Cfcr 



rao e poo 

= / r] sat e br (x)dx + / e br (x)dx (16) 



= A 



%abs 3u ( 

00 p3 

3t _1 e"*d* 

x abs/&e 



bi- 



I e " dt+ 



(17) 



where e^x) is the bremsstrahlung spectrum 



e br (x) = A br exp y~jfj dx ( 18 ) 

and x abs is the absorption frequency u abs in unit of m e c 2 /h. Since A br is smaller by the 
factor exp(x abs /9 e ) in presence of absorption, we take the maximum amplification factor for 
bremsstrahlung as 



v sat = e x ab3 /e e x I e -3 



(19) 



where E n = t n exp(—xt)dt is the exponential integral. The final Comptonized bremsstrahlung 
emission is then 

A Cbr = min(?7o, C*) A br- (20) 
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The energy-weighted mean photon energy for unsaturated Comptonized bremsstrahlung 
is (see e.g. P02) 

4kT x _ J™ ril r xe br (x)dx 



> i-x , , i Vbr@ei 

m e c 2 J r] br e br (x)dx 
and that for saturated Comptonized bremsstrahlung is simply 

T x = T e , 



(21) 



(22) 



since the spectrum approaches Wien spectrum which has Tx = T e . Therefore, the radiation 
temperature of locally radiated bremsstrahlung emission is 



T^ r = -T e mm( Vbr ,4). 



(23) 



The frequency v abs is chosen to be the frequency at which the free-free absorption optical 
depth is equal to 1, 

ra ff (u abs ) = 1, (24) 
where is the absorption coefficient given by Dermer et al. (1991) 



a fA x ) 



o7T- 



ol 2 ott\ 



x 2 el' 2 



1 + (8/7r)V203/ 2 



(25) 



with 



[I + 26 e + 26 2 ) In 



4r7^(l + ZA26 e )e e 



x 



(26) 



+ ( 



26 e )6M 



4^(11.2 + 10.4fl e 2 )fl e 

x 



atf is the fine structure constant, and r e = e 2 /m e c 2 the electron radius (Svensson 1984). 
This expression is valid for i C e and we used (tt /2) 1 / 2 6^ 2 [1 + (S/V) 1 / 2 ^ 2 ] to approximate 
exp(l/0 e )K 2 (l/0 e ). 

Flow at a given position is heated (or cooled) by the inverse Comptonization off elec- 
trons. We use the heating rate (Levich & Sunyaev 1971) 



T c = Aa T c[9 x (r) - 6 e {r)]E x {r)n e {r,&), 



(27) 



where E x {r) is the radiation energy density from equation ([8]) and 9x = kT x /m e c 2 the 
radiation temperature from equation Qj. 



- 9- 



3.2. Equilibrium Temperature and Overheating 

In original CDAF, the temperature of the gas is determined by the balance between the 
viscous heating plus the convective energy transport versus the radiative cooling. However, 
if the radiative heating is dominant in some region of the flow (the condition under which 
this assumption is valid is discussed in §3.4.), the temperature in that region will change to 
reach a new equilibrium. The new equilibrium temperature of the flow will be determined 
by the balance between radiative heating and radiative cooling. 

In CDAF considered in this work, Compton heating versus Comptonized bremsstrahlung 
cooling are the main radiative processes. The thermal equilibrium temperature T eq , then, 
satisfies 

T c (T eq ) = A Cb r(T eq ) (28) 
at given position (r, i?). From equations (I2"0|) and (|2"71) . 9 eq = kT eq /m e c 2 is determined by 

4cE x (r)[6 x (r) -9 eq ] = a f m e c z n e {r,-d)r]cbr{n^r,B eq ){F ei {e eq ) +Fj6 eq )}. (29) 

Since the number density n e (r, $) decrease as $ — > (toward the pole), the derived electron 
temperature increases toward the pole as long as 9x > e - 

For smaller enough the thermal equilibrium temperature can be higher than the virial 
temperature, 

T eq {r,d) > T vir {r). (30) 

The virial temperature is defined as (5/2)kT vir = m p GM/r, where m p is the proton mass. 
Once electrons are heated above the virial temperature, and therefore above the ion tem- 
perature, electrons are likely (via collisions and instabilities) to heat ions to above the virial 
temperature, thereby, creating winds, especially because the dynamical time of CDAF is 
much longer than the free-fall flow. Therefore, we adopt equation (130]) as the condition for 
overheating and producing a wind. 

What equation (I3"U|) means is that some region of CDAF will be radiatively heated from 
the background temperature T e to T eq to achieve the thermal equilibiurm. However, as the 
temperature of the flow approaches T eq , it will become unbound when the condition (|30|) is 
met, and thereby producing the outflow. 

There is also a trivial, yet additional constraint for overheating: the radiation temper- 
ature must be higher than the virial temperature 

T x {r) > T vir (r), (31) 

otherwise the radiation field will cool the flow. 
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3.3. Outflow 

Outflows will extend from the polar axis to the angle $ c at which the equilibrium 
temperature is equal to the virial temperature, 

T eq (r,$ c )=T vir (r). (32) 

Regions with d < i? c are overheated to above the virial temperature. Hence, the shape of 
the outflow is determined by the angle i? c (r) as a function of r. 

One example of the outflow is shown in Figure 3 for the mass accretion rate of m — 10~ 3 , 
the electron temperature of inner region T\ = 10 11 K. The solid curve shows i? c (r) , and the 
outer dotted circle is the outer boundary of the flow. The outflow starts at r ~ 40 and 
the opening angle i9 c reaches maximum of 17° around r ~ 700 and ends with 8.3° at the 
outer boundary. Figure 4 shows the same for rh = 2.0 x 10 -3 and 7\ = 3.0 x 10 10 K with 
i? c = 6.6° at the boundary. The inner isothermal region and the immediate surroundings 
have Tx < T vir , and the overheating does not occur. However, at larger radius Compton 
heating is strong enough to overheat the polar region of the gas and produce outflow. The 
fraction of the sphere covered by the outflowing gas is 27r# 2 /Atx ~ 0.01 for d = 8.3°. 

We have considered combinations of rh and T\. Outflow solutions exist within a limited 
range of these parameters. Too high rh makes the outer boundary so small that the whole 
region is isothermal, i.e., Tx = T e , and the flow is not heated. The outer boundary is smaller 
for higher rh because radiative luminosity is proportional to m 2 while CDAF flows have a 
fixed (assumed) radiation efficiency, which makes the total luminosity roughly proportional 
to rh. The only way to reconcile this is r out being smaller for higher rh (BNQ). In a too 
low rh flow, Comptonization is not strong enough to keep the radiation temperature high. 
Lower energy photons from larger radii dilute the radiation field and flows are not heated. 
In Figure 5, we show region of space in (m,Ti) where outflow solutions are successfully 
produced. Circles denote (m,Ti) for which outflows are found. The size of circle represents 
the opening angle of the outflow at the outer boundary, the majority being < 10°. Crosses 
denote (m,Ti) values for which outflows do not exist. 

As expected, flows with high electron temperature are more likely to develop outflows. 
Outflows are not expected for Ti < m e c 2 /k ~ 6 x 10 9 K: Comptonization is less efficient 
when kT e /m e c 2 < 1. Since the adopted radiation efficiency is e c r] c ~ 10~ 2 , we expect CDAFs 
with outflows to exist in the luminosity range 8 x 10~ 7 < L/Le < 4 x 10~ 5 for T\ = 10 11 K 
and 4 x 10~ 5 < L/L E < 10 -4 for 7\ = 10 10 K. These results do not depend on the mass of 
the black holes. 

So far, we have assumed that the total radiation efficiency e c r] c is constant. In fact, 
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current theoretical understanding of the CDAF solution is not yet secure on this point. 
We now relax that assumption and study an alternative case when e c r] c is assumed to be 
proportional to fa. We adopt simply e c r] c = 10 _2 (m/10 -2 ) so as to make e c r] c = 10~ 2 for 
fa = 10~ 2 . Now for small fa the total luminosity determined from equation (|11|) and the 
total bremsstrahlung emission from equation (fit)]) are both proportional to fa 2 , and the 
outer boundary is now at a roughly constant radius (see dotted lines in Fig. 2). At high 
fa, Comptonization adds additional power of fa to the total emission, and r out depends on 
fa (see dotted lines in Fig. 2). The outflow solutions for this choice of radiation efficiency 
are shown in Figure 6. Qualitatively, they are not much different from the solutions with 
constant e c r] c . Nonetheless, the flow is much less extended due to lower radiation efficiency 
for lower fa and the shape of the outflow funnel in this case is almost straight. The whole 
flow structure is quite self-similar except for the magnitude of the opening angle. Most of 
the solutions still have the opening angle at the outer boundary < 20°. 



3.4. Convection, Heating, and Cooling Timescales 

In the original self-similar CDAF, the viscous heating is balanced by the convective en- 
ergy transport. However, radiative heating can be dominant in some region of the flow. The 
importance of each heating and cooling process is measured by the corresponding timescale. 
The radiative heating timescale is given by the ratio between the internal energy of the 
electron, e e , and the Compton heating rate, 

t H = ^- (33) 
1 c 

The radiative cooling timescale is similarly given by 

to = ( 34 ) 

^Cbr 

The timescale for convective energy transport at a certain two-dimensional position is not 
trivial to evaluate because the two-dimensional velocity profile including the convective mo- 
tion is needed. Although two-dimensional numerical simulations (Stone, Pringle, Sz Begel- 
man 1999; Igumenshchev, Abramowicz, & Narayan 2000; Igumenshchev & Abramowicz 
2000) suggest very limited flow motion near the pole as in self-similar two-dimensional ADAF 
(Narayan & Yi 1995), the two-dimensional flow motion has yet to be expressed in simple 
analytic form. Here, we adopt the height-averaged convective energy flux given by NIA, 

„ c 2 ds , . 
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where a c is the convection coefficient analogous to the usual Shakura & Sunyaev a, c s the 
isothermal sound speed, Qk the Keplerian angular velocity, and s the entropy of the flow. 
In a self-similar solution that is marginally stable to convection, a c = 3a is expected when 
a < 0.05 (NIA). The cooling rate per volume due to this convective energy flux is then 

F c 

Qcv = — . (36) 
r 

Corresponding timescale t cv is now 

t a = (37) 

Since the actual convective motion near the pole is generally much smaller than the average 
convective motion near the disk midplane as is assumed in equation (I35p . the actual timescale 
for the convective cooling could be much longer than t cv evaluated by equation ( 1371) . 

For given dimensionless mass accretion rate rh and the inner electron temperature T%, 
small a (of Shakura-Sunyaev), or equally small a c , means slower radial motion, and therefore, 
slower convective motion. So we searched for the critical value of a below which t cv > t H 
is satisfied within the overheated region. For higher electron temperature of T\ = 10 11 K. 
the value of the critical a ranges from 2.0 x 10~ 3 for rh = 2 x 10~ 3 to 1.2 x 10~ 5 for 
rh = 8 x 10~ 5 , all at a constant radiation efficiency of 10~ 2 . For lower electron temperature 
of T\ = 3 x 10 10 K. the value of the critical a ranges from 2.7 x 10~ 2 for rh = 8 x 10~ 3 to 
3.0 x 10~ 4 for rh = 8 x 10~ 4 . For even lower electron temperature of T\ = 10 10 K at which 
overheated region is very limited in radius, the value of the critical a ranges from 1.6 x 10~ 2 
for rh = 10~ 2 to 2.8 x 10 -3 for rh = 4 x 10~ 3 . Expected values of a for which CDAF exists 
are in the range a < 0.05 (NIA). We, therefore, conclude that the Compton heated outflow 
is expected to be produced in CDAF if the viscosity parameter a is small enough, exact 
value being dependent on the mass accretion rate. However, as already mentioned, the real 
convective motion near the pole is expected to be much smaller than the average convective 
motion in the disk, the Compton heated outflow near the pole may develop even for much 
larger value of a. 

Another timescales we want to check are tjj versus to at the temperature of the back- 
ground flow T e (r,r9) given by equation (jSJ). Only when radiative heating is greater than 
cooling, i.e., tn\T e < tc\r e , the flow is heated from T e to overheated temperature T eq of 
equation (j2"8~|) . The ratio (tn/tc)^ can be expressed as 

t c _ r c (T e ) _ T c {T eq ) r c (T e ) A Cbr (Te q ) _ Tc{T e ) A C br(T eq ) 

tH A C br(T e ) A C br(T eg ) T C (T eg ) A C br{T e ) F c (T eq ) A C br(T e ) 

which is always less than 1 as long as T eq > T e because Y c {T e ) is a monotonically de- 
creasing function of T e and Acbr{T e ) a monotonically increasing function of T e . The ratio 
r 'c{T e q) I Acbr(T eq ) = 1 by the definition of T eq in equation (l28~l) . 
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3.5. Synchrotron Emission 

The synchrotron emission of hot electrons can produce copious soft photons some of 
which are subsequently Compton upscattered. This will lead to an increase in the total 
cooling rate as well as a decrease in the radiation temperature of the emitted radiation. 
Increased cooling rate will reduce the outer boundary of the CDAF whereas decreased ra- 
diation temperature will reduce the Compton heating, or even cool the flow if the radiation 
temperature falls below the electron temperature. 

So we recalculated the condition for overheating in presence of magnetic field for the flow 
parameters considered in Fig. 5. The treatment of Comptonized synchrotron is described in 
Appendix B. We fix the black hole mass to be 1O 8 M . We find that the synchrotron emission 
can significantly affect the high temperature flow whereas the low temperature flow is less 
affected. The overheating and subsequent production of the outflow is still possible for flow 
with T\ = 10 10 K as long as the magnetic field is less than 10% of the equipartition (the gas 
pressure being equal to the magnetic pressure) field. The outer boundary of the flow and the 
opening angle of the outflow are hardly changed if this condition is met. For 7\ = 3 x 10 10 K, 
we find that the magnetic field has to be smaller than 3% of the equipartition field to produce 
the outflow. For the highest temperature 7\ = 10 11 K flow, the production of the outflow 
is possible only when the magnetic field is less than 1% of the equipartition field. So we 
conclude that the Compton preheated outflow is possible in the magnetic CDAFs as long as 
the magnetic field is less than from 10% to 1% of the equipartition field. 



4. Summary and Discussion 

Hot accretion flows like ADAFs have a number of physical characteristics that com- 
pliment the classic low-temperature disk flows. In previous work, we have explored the 
consequences of the high temperature and the two-dimensional density structure of these 
flows, and have found that ADAFs may be able to produce radiatively driven outflows. We 
subsequently have noted that self-similar CDAFs have even more suitable properties for 
producing outflows as compared to ADAFs: steeper poloidal density gradients and higher 
radiation efficiencies. In this paper, we have studied the conditions for self-similar two- 
dimensional CDAFs (NIA; QG) to develop radiatively heated polar outflows. 

1. We have found that CDAFs produce enough luminosity and photon energy to drive 
polar outflows via Compton heating for a reasonable range of mass accretion rate, or, equally, 
luminosity, as long as the magnetic field is less than from 10% to 1% of the equipartition field. 
When the electron temperature saturates around 10 11 K at the inner region, polar outflows 
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are possible for 8 x 1CT 7 < L/L E < 4 x 1CT 5 for radiation efficiency of 10~ 2 , where L E is the 
Eddington luminosity. The luminosity range for which outflow exists is narrower for lower 
electron temperature flows and disappears completely for electron temperature <6x 10 9 K 

2. In most cases, outflows are well collimated along the rotation axis, with an opening 
angle typically in the range < 10°. 

3. If we, instead of taking efficiency as constant, assume that it is proportional to the 
mass accretion rate rh, the solutions are qualitatively the same but are more self similar, i.e., 
opening angle and outer boundary (in Schwarzschild units) depend less strongly on m. 

4. Outflow is more probable for small viscosity parameter a. 

The treatment in this work is not completely satisfactory in the sense that the dynamics 
and the temperature profiles of CDAFs are not self-consistently solved. However, it was not 
our intention to solve fully three-dimensional, self-consistent gloabl CDAFs with proper 
consideration for all gas, radiative, and magnetic processes, which will be eventually needed 
to fully understand CDAFs. Rather, our goal has been to show that even in the framework of 
simple self-similar solutions, radiatively heated polar outflows appear as natural consequences 
of the physical characteristics of CDAFs. 

This work is the result of research activities (Astrophysical Research Center for the 
Structure and Evolution of the Cosmos) supported by Korea Science & Engineering Foun- 
dation. 



A. Radiation Field Inside An Optically Thin Spherical Shell 

Let's consider an optically thin, uniform radiating shell with radius r' and the thickness 
Ar'. The specific intensity at radius r (r < r') in the direction of d (see Figure 7) is given by 
integration of radiative transfer equation along that direction. If we denote the emissivity 
(per unit volume per unit solid angle) of the shell as £ v /4ti, then the specific intensity is 
simply the emissivity times the path length, 

7 " = ^ As = 7^~~~a' ( A1 ) 
47r Air cos 

where is the angle between the ray and the normal of the shell. From the law of sines, 

sin0 = (r/R) sint?. The radiation energy density, E u (r), is given by the integral of I v over 
the solid angle dfl, 

1 f e Ar' 

E v {r) = - / -p— r(Kl (A2) 

C J All COS (J) 
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-Ar' / =dip (A3) 



c Air 



1 - ^ sin 2 



c r 



I 7" I r 1 A- v 

-e.Ar'-lnW-— -. (A4) 



Equation (1A2|) formally diverges when r — > r' whereas the correct value saturates. In the 
limit when r' — r < Ar', the effect of shell's curvature has to be incorporated in As. However, 
since the radiation energy density is calculated by the sum of contributions from discrete 
shells in r (eq. [8]), the divergence is automatically avoided. 



B. Comptonized Synchrotron Emission 



The angle-averaged synchrotron emission by relativistic Maxwellian electrons is given 
by (Pacholczyk 1970) 

2tt e 2 n e v Tll 

where xm = 2v/(3v 9l), u = eB/(2nm e c), and 



e syn {v)dv = ——-^ r l'{x M )dv (Bl) 



r „ , 4.0505 / 0.40 0.5316 \ , , nnnn 1/3 , 

I'M = —n^ 1 + ^71 + -JIT exp(-1.8899xi/ 3 ) (B2) 

nn • \ fTf ' ff ' I 

x Af \ X M X A/ / 

is a fitting formula (Mahadevan, Narayan, & Yi 1996). When absorption is not important, 
the cooling rate due to the optically thin synchrotron emission is obtained by integrating 
the equation (1B1I) 

A% n = 213.6-n e z/ 2 e (B3) 



c 



However, a large fraction of the low energy synchrotron photons are generally absorbed 
by synchrotron self-absorption. The synchrotron emission in the presence of absorption can 
be approximated as 

h-syn = fsynA-syw (B4) 

where 

fsyn = / X M l'(x M )dXM / / X M l'(x M )dx M (B5) 
Jxf} 3 ' JO 

is the fraction of synchrotron emission above the synchrotron self-absorption frequency v a b s 
which satisfies 

1 e 2 cn e {r)r 

T *»" = wi^m 1 M = h (B6) 
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Thus, we only consider the optically thin part of synchrotron emission as the cooling function 
for the gas and as a contribution to the preheating radiation field. 

Locally emitted synchrotron photons are upscattered by inverse Comptonization off hot 
electrons. Here, we adopt a simple estimate of Comptonized synchrotron, which is reasonable 
in physical conditions considered in this paper, 

Acs = VsynfsynKym ( B7 ) 

where 

7] syn = min(^ ,<n) (B8) 

is the Comptonized energy enhancement factor (see e.g., Dermer et al. 1991) for synchrotron 
emission. The enhancement factor for fully saturated Comptonized synchrotron in the pres- 
ence of absorption is 

Labs ^f~%Ml'{ x M)dXM 

rill = % T „ w > (B9) 



f^bs x M I'{x M )dxM 



A/ 

(see §3.1 for bremsstrahlung case). The synchrotron absorption frequency is compared with 
that for free-free absorption, and the larger of the two is chosen. 

The radiation temperature of locally emitted Comptonized synchrotron is similarly 

T% s = min^o^^ 2 <x M >e v ,T e ), (BIO) 

where 

Labs x M I\xM)dxM 

< X M >e = -j£ „, w (Bll) 

Lab* x M i {x M )dx M 

JX M 

is the energy-weighted mean photon energy in xm unit. This treatment of Comptonized 
synchrotron radiation is almost the same as that adopted in Park & Ostriker (2001), except 
that the energy enhancement factor and the treatment of saturated Comptonization are 
slightly improved. 

In the presence of both bremsstrahlung and synchrotron emission, A is replaced by 
A C br + Acs in equation (|S]) and T X A by T^ br A C br + T^ s A C s m equation 
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Fig. 1. — Typical log profiles of luminosity L(r) (in units of Le, left ordinate), of the radiation 
energy density Ex{r) (in arbitrary units, left ordinate), and of the radiation temperature 
T x {r) (in K, right ordinate). 
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Fig. 2. — The outer boundary radius r out as a function of m. Upper solid line is for T\ = 
3 x 10 10 K and lower one for 7\ = 10 11 K, both for constant radiation efficiency e c i] c = 0.01. 
Dotted lines represent the same r out for varying radiation efficiency e c n c = 10~ 2 (m/10" 2 ): 
upper one for T\ = 3 x 10 10 K and lower one for 7\ = 10 11 K. 
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Fig. 3. — The region inside the solid curve toward the pole is overheated above the virial 
temperature due to Compton heating. Outer dotted circle shows the outer boundary of 
CDAF for given total luminosity. This figure is for m = 10~ 3 and T\ = 10 11 K. The opening 
angle of the outflow at the outer boundary is only 8°. 
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Fig. 4. — Same as Figure 3, but for m = 2 x 10 3 and 7\ = 3 x 10 10 K. The opening angle 
at the outer boundary is 8°. 
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Fig. 5. — The parameter space of (m, T\) under constant radiation efficiency for which 
outflow exists is denoted as circles. The size of circle represents the opening angle of the 
outflow at the outer boundary. No outflow is found at (m, Ti) denoted by crosses. 
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Fig. 6. — Same as Figure 5, but now the radiation efficiency is proportional to the mass 
accretion rate. 




Fig. 7. — The diagram for calculation of specific intensity, I v , inside a radiating spherical 
shell. 



